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"(N ' ABSTRACT 
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The temporal and spectral properties of terrestrial gamma-ray flashes (TGFs) are 
\ studied. The delay of low energy photons relative to high energy ones in the 7-ray 

Li^ \ variations of the TGFs with high signal to noise ratio has been revealed by an approach 

I/-) ' of correlation analysis in the time domain on different time scales. The temporal 

structures of the TGFs in high energy band are usually narrower than that in low 
energy band. The spectral hardness has a general trend of decreasing with time in a 
flash. The observed temporal and spectral characteristics give constraints and valuable 
hints on the flash production mechanism. 
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1. Introduction 

Intense 7-ray bursts of atmospheric origin were detected by the Burst and Transient Source 
Experiment (BATSE) detectors, located on the Compton Gamma Ray Observatory (CGRO) 
( Fishman et ai, 1994| ). All the observed TGFs were of short duration (a few milliseconds) and 



CO 



thought to be caused by high-altitude discharges produced by runaway air breakdown. Several 
theoretical calculations have been done to explain the origins of TGFs under the runaway air 
breakdown mechanism ( \Chang and Price, 1995 ; Taranenko and Roussel-Dupre, 1996; Lehtinen 



had extremely hard energy spectra consistent with bremsstrahlum radiation from energetic (MeV) 
electrons. Since they were spatially correlated with regions of high thunderstorm activity and 
^ ■ in two cases have been correlated with individual lightning flashes ( Inan et ai, 1996| ), they are 



et al, 1996 ; Lehtinen et al, 1997 ). The calculated 7-ray fluxes and duration of the bursts are 



comparable to those measured by the BATSE detector on the CGRO satellite. 

As the runaway air breakdown is a fundamental new process in plasma, it is interesting to 
understanding more the physics of the process through studying TGFs, if they were truly caused 
by such a mechanism. Besides the total fluxes, durations and time-averaged spectra, the energy 
dependence of temporal profiles and spectral evolution of TGFs can give more information about 
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their production process. For this purpose we analyze the time-tagged event (TTE) data for 
TGFs observed by BATSE. For 4 energy channels with boundaries 25-55 keV, 55-110 keV, 110-320 
keV, and > 320 keV, the TTE data contains the arrival time (2/xs resolution) of each photon 
from a triggered event up to ~ 3.3 x 10 4 photons. Of 47 TGFs observed on CGRO in ~ 4 years 
and published by NASA in their details, 15 bursts with high signal to noise ratio are used in our 
analysis, which are selected by the criterion of F m /TgQ > 4 ph/cm _2 /0.1ms/ms, where F m is the 
peak rates of the counting series with time step 0.1 ms, T90 (ms) the burst duration which contains 
%90 burst photons. 



2. Energy Dependence of Temporal Profiles 

To study the energy dependence of temporal profiles, which is important to diagnose the 
emission mechanism, we construct two counting series /i(t) and f2(t) for each burst for the energy 
band 25-110 keV and > 110 keV, respectively. The average pulse widths in the low and high 
energy bands and the relative time delays between the two bands can be calculated by a modified 
correlation analysis technique used in studying X-ray rapid variability of the black hole binary 
Cyg X-l ( gj et al, 1999| ; \Li, 2001Q . 



The cross-correlation function of two time series /1 and ji at time lag r is defined as 

CCF(r) = J2M^t + r)v2(iAt)/a( Vl )a(v 2 ) (1) 

i 

where v{t) = f{t) — f, f{t) is the number of photons in the time interval (t, t + At), At is the time 
step. If the function CCF(r)/CCF(0) has maximum at r = A, the time lag of the energy band 1 
relative to the band 2 at the time scale Deltat is then defined as A. Monte Carlo simulations have 
been done and the results show that with this technique we can measure the relative time delay 
between two bands over a wide range of time scale At with high time resolution. At large scales 
one usually can get enough signal photons in a time bin and desirable correlation values from finite 
time bins. And at small scales the effect of serious Poisson fluctuation of signal counts in a time 
bin can be compensated by the large amounts of time bin and accurate correlation values can also 
be derived with Eq. (1). A distribution of time lag vs. time scale can reflect the character of the 
physical process to produce the delay better than a single value of lag at only one time scale. A 
physical process usually occurs in a range of time scale, the spectral delay caused by the process 
should appear at different time scales, smoothly distributed in the range. On the other hand 
apparent delays from statistical fluctuation will fluctuate between positive and negative values. 

For each studied burst and m = 25 different values of time step At which are logarithmically 
uniformly placed in the region of 10 -5 — 10 -3 s, we calculate the time lags A of f± relative to 
fi- All obtained time lags Aj (i = l,...,m) of each burst with high signal to noise ratio are 
always positive. We average each 5 successive Aj and show the distribution of average time 
lag vs. time scale in Fig. 1. The global average A = Y^=\^il m an d the standard deviation 
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cr(A) = y X)(Aj — A) 2 /(m — 1) for each selected TGF is listed in Table 1. The total counting 
profile and profiles in 4 energy bands of a TGF with BATSE trigger number 2955 are plotted in 
Figure 2, where the delay of lower energy photons relative to higher energy ones is apparent. 

The width W\ of a temporal profile in a band I can be defined as the FWHM of the 
autocorrelation function 

ACF(r) = £ Vl (iAt + T>,(iAt)/c7 2 («i) (2) 

i 

The widths W\ and W% of studied TGFs in the low and high energy bands are calculated and 
their ratios W1/W2 are presented in Table 1, from which we can see that the lower energy pulses 
are wider than higher energy ones for most studied bursts. 



3. Spectral Evolution 

For two counts c\ in 25-110 keV band and C2 in > 110 keV band recorded in the same time 
segment, we defined the hardness ratio as 

h = ^^ (3) 
c 2 + ci 

To show the spectral evolution, a TGF's pulse is divided into 5 time segments, each segment 
contains approximately the same number of photons and the hardness ratio is calculated by Eq. 
(3). Figure 3 shows the spectral hardness variation vs. time (indicating by the segment number) 
with four panels grouped by the shape of temporal evolution of hardness. Although the shapes are 
various, a general trend of hardness that decreases with time appears in Fig. 3, which confirms 
the claim by Nemiroff et al. [1997] with earlier TGF data. 



4. Discussion 

For all the bursts with high signal to noise ratio we find that in comparison with high energy 
band of > 110 keV, 7-ray variations in the low energy band of 25 - 110 keV are always late in the 
order of ~ 100 [as, pulses are usually wide, and the energy spectra have a general trend to softening 
with time. The above features of energy dependence of time profiles and spectral revolution 
observed in TGFs support models of runaway breakdown as they are naturally expected for 7-rays 
produced by explosive discharges in plasmas. Monte Carlo calculations show that large-voltage 
and high-temperature pinch plasma columns can generate observed 7-ray flashes with energy 
spectra and spectral evolution characteristics, including the magnitude of soft lags, consistent with 



those observed in 7-ray bursts ( Li and Wu, 1997 ; Li, 199 



Alternative mechanisms have been proposed to explain the origin of TGFs, e.g. Fargion [2001] 
claimed that TGFs originated from 10 15 — 10 eV neutrinos of astrophysical nature induced air 
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showers and take TGFs as first evidences for discovering ultra high energy (UHE) neutrinos. Our 
Monte-Carlo simulations show that up-going air showers are hard to produce the observed time 
lags between hard and soft energy photons in TGFs. In our simulations a geometrical model of 
the earth, atmosphere and BATSE detector is set. High energy 7-rays initiated by UHE neutrinos 
are produced from the earth's surface, following with the up-going air showers. We tracked the 
secondary photons and electrons and record the arriving time and energy of all 7-rays reaching 
the detector, with deposition energy in the BATSE's region. Simulations were made with the 
Monte-Carlo package GEANT3 from CERN. In more than one thousand simulated bursts no time 
lag between hard and soft photons larger than 10 //s has been found. 

Quantitative comparison between observed features and expectations from models are needed 
to further confirm the mechanism and to study physics of the production process. 

We thank Prof. Y. Muraki and the reviewers for helpful comments and suggestions. This 
study was supported by the Special Funds for Major State Basic Research Projects and by the 
National Natural Science Foundation of China and made use of data obtained through HEASARC 
Online Service, provided by NASA/GSFC. 
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Table 1 Selected TGFs 



Trig No. 


T 90 (ms) 


F m (10~ 3 ph/cm 2 /0.1ms) 


Fm/Tgo 


A±cr(A) (fis) 


W 1 /W 2 


106 


1.3 


6.74 


5.19 


102 ± 28 


1.19 


1433 


1.5 


7.15 


4.77 


71 ± 31 


1.64 


2144 


0.7 


9.62 


13.74 


115 ± 40 


1.03 


2348 


1.1 


7.15 


6.50 


88 ± 29 


1.04 


2370 


0.7 


8.14 


11.63 


151 ± 29 


1.16 


2465 


1.0 


7.89 


7.89 


157 ± 25 


1.38 


2754 


0.9 


7.65 


8.50 


120 ± 29 


2.10 


2808 


1.0 


5.67 


5.67 


105 ± 45 


0.88 


2955 


0.7 


7.89 


11.28 


125 ± 21 


1.00 


3377 


1.0 


10.11 


10.11 


373 ± 31 


1.47 


5577 


0.6 


5.18 


8.63 


81 ± 33 


0.91 


5587 


1.2 


6.41 


5.34 


163 ± 46 


1.16 


5598 


1.1 


6.17 


5.61 


117 ± 20 


0.94 


5665 


1.4 


7.65 


5.46 


265 ± 43 


1.54 


6773 


0.9 


3.82 


4.25 


278 ± 61 


0.74 
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Fig. 1. — Soft 7-ray lag vs. time scale of TGFs. The quantity beside an arrow is the BATSE trigger 
number of the indicated burst. 




Fig. 2. — Counting rate profiles of a TGF with trigger number 2955 in different energy bands. 
From top to bottom is the profile in 25-60 keV, 60-110 keV, 110-325 keV, > 325 keV and the total 
BATSE band, respectively. 
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Fig. 3. — Spectral evolution of TGFs. The vertical axis is the hardness ratio defined by Eq. (3). 
The number marked along the abscissa is the segment number of a burst profile 



